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Abstract. 

The DAMA/LIBRA set-up (about 250 kg highly radiopure Nal(Tl) sensitive mass) is running at 
the Gran Sasso National Laboratory of the I.N. F.N. . The first DAMA/LIBRA results confirm the 
evidence for the presence of a Dark Matter particle component in the galactic halo, as pointed out 
by the former DAMA/Nal set-up; cumulatively the data support such evidence at 8.2 o C.L. and 
satisfy all the many peculiarities of the Dark Matter annual modulation signature. The main aspects 
and prospects of this model independent experimental approach will be outlined. 
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INTRODUCTION 

DAMA is an observatory for rare processes located deep underground at the Gran 
Sasso National Laboratory of the I.N.F.N.. It is based on the development and use of 
low background scintillators. The main experimental set-ups are: i) DAMA/Nal (~ 
100 kg of highly radiopure Nal(Tl)) which completed its data taking on July 2002 
CO, 0, 1, H SlEE llEEE Q; m DAMA/LXe (~ 6.5 kg liquid Kr-free Xenon 
enriched either in 129 Xe or in 136 Xe) lfl3|l : iii) DAM A/R&D, a facility dedicated to tests 
on prototypes and to experiments for investigations on rare events | j 14T | : iv) DAMA/Ge, 
dedicated to sample measurements and to specific measurements on rare events lfl5ll : 
y) the second generation DAMA/LIBRA set-up (~ 250 kg highly radiopure Nal(Tl)) 
Profiting of the low background features of these set-ups, many rare 
processes are studied. 

In particular, DAMA/LIBRA is investigating the presence of Dark Matter (DM) 
particles in the galactic halo by exploiting the model independent DM annual modulation 
signature. This signature - originally suggested in the middle of 80's [19] - exploits the 
effect of the Earth revolution around the Sun on the number of events induced by the 
Dark Matter particles in a suitable low-background set-up placed deep underground. In 
fact, as a consequence of its annual revolution, the Earth should be crossed by a larger 
flux of Dark Matter particles around ~ 2 June (when its rotational velocity is summed 



to the one of the solar system with respect to the Galaxy) and by a smaller one around 
~ 2 December (when the two velocities are subtracted). This offers an efficient model 
independent signature, able to test a large number of DM candidates, a large interval of 
cross sections and of halo densities. It should be stressed that the DM annual modulation 
is not - as often naively said - a "seasonal" variation and it is not a "winter- summer" 
effect; in fact, the DM annual modulation is related to the Earth velocity in the galactic 
frame and its phase (roughly 2 nd June) is well different than those of physical quantities 
(such as temperature of atmosphere, pressure, other meteorological parameters, cosmic 
rays flux, ...) correlated with seasons instead. 

The DM annual modulation signature is very distinctive since the corresponding 
signal must simultaneously satisfy all the following requirements: the rate must contain 
a component modulated according to a cosine function (1) with one year period (2) and 
a phase that peaks roughly around ~ 2 nd June (3); this modulation must only be found 
in a well-defined low energy range, where DM particle induced events can be present 
(4); it must apply only to those events in which just one detector of many actually "fires" 
(single-hit events), since the DM particle multi-interaction probability is negligible (5); 
the modulation amplitude in the region of maximal sensitivity must be <7% for usually 
adopted halo distributions (6), but it can be larger in case of some possible scenarios 
such as e.g. those in refs. [|20L \2l \ j. Only systematic effects or side reactions able to fulfil 
these requirements and to account for the whole observed modulation amplitude could 
mimic this signature; thus, no other effect investigated so far in the field of rare processes 
offers a so stringent and unambiguous signature. 

At present status of technology it is the only model independent signature available in 
direct dark matter investigation that can be effectively exploited. 

It is worth noting that the corollary questions related to the exact nature of the DM 
particle(s) (detected by means of the DM annual modulation signature) and to the as- 
trophysical, nuclear and particle Physics scenarios require subsequent model dependent 
corollary analyses, as those performed e.g. in refs. JIB S 0, S B OS H • On the other 
hand, one should stress that it does not exist any approach in direct and indirect DM 
searches which can offer information on the nature of the candidate in a model inde- 
pendent way, that is without assuming some astrophysical, nuclear and particle Physics 
scenarios. 

In the following, we will just briefly summarize the first results on the Dark Matter 
particle investigation obtained by DAMA/LIBRA, exploiting over four annual cycles 
the model independent DM annual modulation signature (exposure of 0.53 tonxyr). 
The results have also been combined together with the previous data collected over 7 
annual cycles by DAMA/Nal (0.29 tonxyr). Thus, the whole available data correspond 
to 11 annual cycles for a total exposure of 0.82 tonxyr, which is orders of magnitude 
larger than the exposures typically collected in the field. 



DAMA/LIBRA RESULTS 



The DAMA/Nal set up and its performances are reported in ref.[Ql |3j3, 01, while 
the DAMA/LIBRA set-up and its performances are described in ref. 111611 . Here we 
just summarized the main features: i) the sensitive part of the set-up is made of 25 



highly radiopure Nal(Tl) crystal scintillators placed in a 5-rows by 5-columns matrix; 
ii) the detectors' responses range from 5.5 to 7.5 photoelectrons/keV; iii) the hardware 
threshold of each PMT is at single photoelectron (each detector is equipped with two 
low background photomultipliers working in coincidence); iv) energy calibration with 
X-rays/y sources are regularly carried out down to few keV in the same conditions as 
the production runs; v) the software energy threshold of the experiment is 2 keV. 
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FIGURE 1. Experimental model-independent residual rate of the single-hit scintillation events, mea- 
sured by DAMA/Nal and DAMA/LIBRA in the (2 - 6) keV energy interval as a function of the time. The 
zero of the time scale is January I st of the first year of data taking of DAMA/Nal. The experimental points 
present the errors as vertical bars and the associated time bin width as horizontal bars. The superimposed 
curve is the cosinusoidal functions behaviors A cos co(t — to) with a period T = ^2 = 1 yr, with a phase 
fo = 152.5 day (June 2 nd ) and with modulation amplitude, A, equal to the central value obtained by best 
fit over the whole data, that is: (0.0129 ± 0.0016) cpd/kg/keV. The dashed vertical lines correspond to the 
maximum of the signal (June 2 nd ), while the dotted vertical lines correspond to the minimum. The total 
exposure is 0.82 ton x yr. For details and for other energy intervals see 11711 . 

Several analyses on the model-independent DM annual modulation signature have 
been performed (see ref. [17] and references therein); here just few arguments are men- 
tioned. In particular, Fig. Q] shows the time behaviour of the experimental residual rates 
of the single-hit events collected by DAMA/Nal and by DAMA/LIBRA in the (2-6) keV 
energy interval 11711 . The superimposed curve in Fig.[TJrepresents the cosinusoidal func- 
tions behavior Acos co(t — to) with a period T = = 1 yr and with a phase to = 152.5 
day (June 2 nd ), while the modulation amplitude, A, has been obtained by best fit over 
the DAMA/Nal and DAMA/LIBRA data. When the period and the phase parameters 
are also released in the fit, values well compatible with those expected for a DM particle 
induced effect are obtained lfl7tl: T = (0.998 ± 0.003) yr and t = (144 ± 8) day in the 
cumulative (2-6) keV energy interval. Summarizing, the analysis of the single-hit resid- 
ual rate favours the presence of a modulated cosine-like behaviour with proper features 
at8.2aC.L. |Q3]. 

The same data of FigQ]have also been investigated by a Fourier analysis, obtaining 
a clear peak corresponding to a period of 1 year QJ/ZU ; this analysis in other energy 
region shows instead only aliasing peaks. Moreover, while in the (2-6) keV single- 
hit residuals a clear modulation is present, it is absent at energies just above DJ/7j]. In 
particular, in order to verify absence of annual modulation in other energy regions and, 
thus, to also verify the absence of any significant background modulation, the energy 
distribution measured during the data taking periods in energy regions not of interest 



for DM detection has also been investigated. In fact, the background in the lowest 
energy region is essentially due to "Compton" electrons, X-rays and/or Auger electrons, 
muon induced events, etc., which are strictly correlated with the events in the higher 
energy part of the spectrum; thus, if a modulation detected in the lowest energy region 
would be due to a modulation of the background (rather than to a signal), an equal or 
larger modulation in the higher energy regions should be present. The data analyses 
have allowed to exclude the presence of a background modulation in the whole energy 
spectrum at a level much lower than the effect found in the lowest energy region for the 
single-hit events II 1711 . 

A further relevant investigation has been done by applying the same hardware and 
software procedures, used to acquire and to analyse the single-hit residual rate, to the 
multiple-hits one. In fact, since the probability that a DM particle interacts in more than 
one detector is negligible, a DM signal can be present just in the single-hit residual rate. 
Thus, this allows the test of the background behaviour in the same energy interval of the 
observed positive effect. In particular, Fig. [2] shows the residual rates of the single-hit 
events measured over the four DAMA/LIBRA annual cycles, as collected in a single 
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FIGURE 2. Experimental residual rates over the four DAMA/LIBRA annual cycles for single-hit events 
(open circles) (class of events to which DM events belong) and for multiple-hits events (filled triangles) 
(class of events to which DM events do not belong), in the energy interval (2 - 6) keV. They have been 
obtained by considering for each class of events the data as collected in a single annual cycle and by 
using in both cases the same identical hardware and the same identical software procedures. The initial 
time of the scale is taken on August 1 th . The experimental points present the errors as vertical bars and 
the associated time bin width as horizontal bars. See ref. 11711 . Analogous results were obtained for the 
DAMA/NaldataJl. 

annual cycle, together with the residual rates of the multiple-hits events, in the same 
considered energy interval. A clear modulation is present in the single-hit events, while 
the fitted modulation amplitudes for the multiple-hits residual rate are well compatible 
with zero 0/7] . Similar results were previously obtained also for the DAMA/Nal case ^ . 
Thus, again evidence of annual modulation with proper features, as required by the DM 
annual modulation signature, is present in the single-hit residuals (events class to which 
the DM particle induced events belong), while it is absent in the multiple-hits residual 
rate (event class to which only background events belong). Since the same identical 
hardware and the same identical software procedures have been used to analyse the two 
classes of events, the obtained result offers an additional strong support for the presence 
of a DM particle component in the galactic halo further excluding any side effect either 
from hardware or from software procedures or from background. 
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The annual modulation present at low energy has also been analyzed by depicting 
the differential modulation amplitudes, S m *, as a function of the energy (the k index 
identifies the energy interval); the S m * is the modulation amplitude of the modulated 
part of the signal obtained by maximum likelihood method over the data, considering 
T = 1 yr and to = 152.5 day. The S m ^ values are reported as function of the energy 
in ref. UM- It can be inferred that a positive signal is present in the (2-6) keV energy 
interval, while S OT * values compatible with zero are present just above; in particular, the 
S OT jfe values in the (6-20) keV energy interval have random fluctuations around zero with 
X 2 equal to 24.4 for 28 degrees of freedom. It has been also verified that the measured 
modulation amplitudes are statistically well distributed in all the crystals, in all the 
annual cycles and energy bins; these and other discussions can be found in ref. Ill7n . It 
is also interesting the results of the analysis performed by releasing the assumption of a 
phase to = 152.5 day in the procedure of maximum likelihood to evaluate the modulation 
amplitudes from the data of the seven annual cycles of DAMA/Nal and the four annual 
cycles of DAMA/LIBRA. In this case alternatively the signal has been written as: 
So,k + S m ,kcosco(t -t ) + Z m ^ sin co(t -to) = So,* + Y m ,k cos CD (*-**), where So,* is the 
constant part of the signal in k-th energy interval. Obviously, for signals induced by DM 
particles one would expect: i) Z m £ ~ (because of the orthogonality between the cosine 
and the sine functions); ii) S w * ~ F w *; iii) t* — to = 152.5 day. In fact, these conditions 
hold for most of the dark halo models; however, it is worth noting that slight differences 
in the phase can be expected in case of possible contributions from non-thermalized 
DM components, such as e.g. the SagDEG stream [0] and the caustics ||2~2ll . The 2a 
contours in the plane (S m ,Z m ) for the (2-6) keV and (6-14) keV energy intervals and 
those in the plane (Y m ,t*) are reported in ref. |17]. The best fit values for the (2-6) keV 
energy interval are (la errors): S m = (0.0122 ±0.0016) cpd/kg/keV; Z m = -(0.0019 ± 
0.0017) cpd/kg/keV; Y m = (0.0123 ±0.0016) cpd/kg/keV; t* = (144.0±7.5) day; while 
for the (6-14) keV energy interval are: S m = (0.0005 ± 0.0010) cpd/kg/keV; Z m = 
(0.0011 ±0.0012) cpd/kg/keV; Y m = (0.0012 ±0.0011) cpd/kg/keV and t* obviously 
not determined. These results confirm those achieved by other kinds of analyses. In 
particular, a modulation amplitude is present in the lower energy intervals and the period 
and the phase agree with those expected for DM induced signals. For more detailed 
discussions see ref. ft 1 VII 

Both the data of DAMA/LIBRA and of DAMA/Nal fulfil all the requirements of the 
DM annual modulation signature. 

As previously done for DAMA/Nal [0, HI], careful investigations on absence of any 
significant systematics or side reaction effect in DAMA/LIBRA have been quantitatively 
carried out and reported in details in ref. 1170. In order to continuously monitor the 
running conditions, several pieces of information are acquired with the production data 
and quantitatively analyzed. No modulation has been found in any possible source of 
systematics or side reactions for DAMA/LIBRA as well; thus, cautious upper limits 
(90% C.L.) on the possible contributions to the DAMA/LIBRA measured modulation 
amplitude have been estimated QJ/7J]- No systematics or side reactions able to mimic 
the signature (that is, able to account for the measured modulation amplitude and 
simultaneously satisfy all the requirements of the signature) has been found or suggested 
by anyone over more than a decade. For detailed quantitative discussions on all the 



related topics and for results see ref. [ITvh and refs. therein. Just as an example we 
recall here the case of muons, whose flux has been reported by the MACRO experiment 
to have a 2% modulation with phase around mid-July [23J]. In particular, it has been 
shown that not only this effect would give rise in the DAMA set-ups to a quantitatively 
negligible contribution lfl7l . 3, @], but several of the six requirements necessary to 
mimic the annual modulation signature - namely e.g. the conditions of presence of 
modulation just in the single-hit event rate at low energy and of the phase value - would 
also fail. Moreover, even the pessimistic assumption of whatever hypothetical (even 
exotic) possible cosmogenic product - whose decay or de-excitation or whatever else 
might produce: i) only events at low energy; ii) only single-hit events; iii) no sizeable 
effect in the multiple-hits counting rate - cannot give rise to any side process able to 
mimic the investigated DM signature. In fact, not only this latter hypothetical process 
would be quantitatively negligible fl7|l . but in addition its phase - as it can be easily 
derived - would be (much) larger than July 15th, and therefore well different from the 
one measured by the DAMA experiments and expected by the DM annual modulation 
signature (~ June 2nd). Recently, a LVD analysis MM has been reported for the muon 
flux relatively to the period 2001-2008, which partially overlaps the DAMA/Nal running 
periods and completely those of DAMA/LIBRA. A value of ~ 185 days has been 
measured by LVD in this period for the muon phase to be compared with (144± 8) 



days [17] which is the measured phase by the DAMA/Nal and DAMA/LIBRA for the 
low energy peculiar single-hit rate modulation. Thus, the latter one is > 5<7 far from the 
muon modulation phase measured at LNGS by the large surface apparata MACRO and 
LVD. In conclusion, any possible effect from muons can be safely excluded on the basis 
of all the given quantitative facts (and just one of them is enough). 

Let us here address in details some recent confusion about the presence of potassium 
in the detectors and its hypotethical role. As first we remark that the only potassium 
isotope contributing to the background is the radioactive 40 K (natural abundance 1 . 17 x 
10~ 4 and half life 1.248 x 10 9 yr), which is at ppt (10~ 12 g/g) level in the detectors 



111611 ; it is worth noting that it appears difficult to do better in Nal(Tl), considering e.g. 
the chemical affinity of Na and K. Recently, it was claimed [25] that some role might 
be played by 40 K; this possibility is evidently discarded by the data, by the published 
analyses and also by simple considerations. Let us summarize that: 

• although the peak around 3 keV in the cumulative energy spectrum (see Fig. 1 of 
ref. fl7|l and the discussion in ref. ||26|l ) can be partially ascribed to 40 K decay, 
there is not evidence for any 3 keV peak in the S ra distribution (see Fig. 9 of ref. 
fl7|l ). At the present level of sensitivity the S m behaviour is compatible within 
the uncertainties both with a monotonic behaviour and with a kind of structure, 
as expected for many Dark Matter candidates and also for WIMPs scenarios; see 
e.g. the Appendix of ref. (|l7l] : 

• 40 K decay cannot give any modulation at all, as well known, unless evoking new 
exotic physics (see later); 

• no modulation has been observed in other energy regions where 40 K decays also 
contribute Ill7ll ; 

• no modulation has been observed in multiple-hits events (events where more than 
one detector fires) in the same energy region where DAMA observes the peculiar 



modulation of the single-hit events (events where just one detector fires). In fact, 
40 K can also give rise to double events in two adjacent detectors when: i) 40 K 
decays in a detector, A, by EC of K shell to the 1461 keV level of 40 Ar; ii) the 
1461 keV 7 escapes from the A detector and hits an adjacent one causing a double 
coincidence. The 3.2 keV X-rays/Auger electrons from K shell of 40 Ar are fully 
contained in the A detector with efficiency ~ 1, giving rise in A to a 3.2 keV peak. 
These double coincidence events, shown in Fig. [3j and also possible multi-site 
events due to Compton scatterings, are multiple-hits events, and are not modulated, 
as fore-mentioned (see also ref. UJ\\ ); 

no modulation is present just in the double coincidence events (3.2 keV - 1461 
keV) due to 40 K decay. In fact, their residuals Q - normalized to the error - as 
function of the time are shown in Fig. ^[bottom left) throughout the four yearly 
cycles of DAMA/LIBRA [17]. The two superimposed curves represent the best-fit 
modulation behaviours Acosco(t — to) and Asinco(t — to) with 1 year period (CO — 
lit /\y ear) and phase to = June 2nd; the best fit values of the modulation amplitudes 
are: A = — (0.10± 0. 12) and A = -(0.02 ± 0. 12), respectively. These latter values 
are well compatible with absence of any annual modulation. Moreover, the data 
show a random fluctuation around zero (% 2 jdof = 1 .08) and their distribution (see 
Fig. [3]) is well represented by a gaussian with r.m.s. ~ 1 . Hence, in conclusion, there 
is absence of any annual modulation in the double coincidence events due to 40 K 
decay; 

fixing the phase of cosine to the perihelion (~ Jan 3rd) in the best fit of the data 
points of Fig. ^bottom left), the modulation amplitude is compatible with zero: 
A = (0.10 ±0.12). Let us note that the effect in ref. [Hzl] - where a 0.3% yearly 
modulation of nuclear decays of two nuclides has been reported H with a phase 
roughly equal to perihelion - is well below this sensitivity, and this effect is not 
able to account for the DAMA signal because of the well-different phase, of the 
marginal modulation amplitude and of several other arguments given in this list, 
the analysis of 40 K double coincidences also rules out at level more than 10c C.L. 
the modulation amplitudes expected when assuming that the DAMA effect might 
be due to the two hypothetical cases of: i) 40 K "exotic" modulation decay (also see 
above); ii) spill-out from double to single events and viceversa along the time; on 
the other hand these exotic hypotetical arguments were already excluded also by 
several other arguments given in this list. 

the behaviour of the overall efficiency during the whole data taking period is highly 
stable. A quantitative investigation of the role of the efficiency as possible system- 
atics leads to a cautious upper limit less than 1% of the modulation amplitude, as 
reported e.g. in II 1711 : 



They are calculated averaging - over the available pairs of adjacent detectors in DAMA/LIBRA - the 
counting rate of the double coincidence events (3.2 keV - 146 1 keV) due to 40 K decay in every pair, once 
subtracted its mean value. 

2 Let us remark that in ref. [28] this effect has already been confuted by searching for modifications to the 
exponential radioactive decay law with the Cassini spacecraft. 
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FIGURE 3. Top: scatter plot of double coincidence events between a DAMA/LIBRA detector, A, 
(low energy region) and an adjacent detector (higher energy region). The threshold of each PMT is at 
single photoelectron level. For comparison, the software energy threshold used in the data analyses of the 
single-hit events for Dark Matter particle investigation: 2 keV, is shown as continuous line. The double 
coincidence events due to 40 K decay (3.2 keV - 1461 keV) are well identified. For details see ref. lfl6tl . 
Bottom left: residuals - normalized to the error - of the double coincidence events (3.2 keV - 1461 keV) 
due to 40 K decay as function of the time throughout the four yearly cycles of DAMA/LIBRA IU7I1 . The two 
superimposed curves represent the best-fit modulation behaviours Acosco(t — to) and As in CO (t — to) with 1 
year period (co = 2n/\year) and phase fo = June 2nd; the best fit values of the modulation amplitudes are: 
A = —(0.10 ±0.12) and A = —(0.02 ±0.12), respectively. Both values are well compatible with absence 
of any annual modulation in the double coincidence events due to 40 K decay. Moreover, the data show 
a random fluctuation around zero (% 2 /dof = 1.08) and their distribution (bottom right panel) is well 
represented by a gaussian with r.m.s. ~ 1. 



the annual modulation signal - observed by DAMA - is present both in the outer 
and in the inner detectors [17] . Hence, there is no dependence on the veto capability, 
that is different - by geometrical reasons - among the outer and the inner detectors. 
In particular, if the 40 K decay would hypothetically play some role in the annual 
modulation of the low-energy single-hit events, the effect would be larger in the 
outer detectors where the 1461 keV y's accompanying the 3.2 keV X-rays/Auger 
electrons have lower probability to be detected by closer detectors; 
the annual modulation signal - observed by DAMA - is equally distributed over all 
the detectors II 1711 : see also the previous considerations. 



Thus, for all the above reasons (and just one of them is enough), no role can be played 
by 40 K. 

In conclusion, DAMA/LIBRA has confirmed the presence of an annual modulation 
satisfying all the requirements of the DM annual modulation signature, as previously 
pointed out by DAMA/Nal; in particular, the evidence for the presence of DM particles 
in the galactic halo is cumulatively supported at 8.2 o C.L.. 

It is worth noting that no other experiment exists, whose result can be directly com- 
pared in a model-independent way with those by DAMA/Nal and DAMA/LIBRA. 
Moreover, concerning those activities claiming for some model dependent exclusion un- 
der some set of largely arbitrary assumptions (see for example discussions in pl l29l . [30|l '). 
some important critical points exist in some of their experimental aspects (energy thresh- 
old, energy scale, multiple selection procedures, disuniformity of the detectors response, 
absence of suitable periodical calibrations in the same running conditions and in the 
claimed low energy region, stabilities, etc.). 

Finally, as regards the indirect detection searches, let us note that no direct model- 
independent comparison can be performed between the results obtained in direct and 
indirect activities, since it does not exist a biuni vocal correspondence between the 
observables in the two kinds of experiments. Anyhow, if possible excesses in the positron 
to electron flux ratio and in the 7 rays flux with respect to an assumed simulation of the 
hypothesized background contribution, which is expected from standard sources, might 
be interpreted in terms of Dark Matter (but huge and still unjustified boost factor and 
new interaction types are required), this would also be not in conflict with the effect 
observed by DAMA experiments. 



UPGRADES AND PROSPECTS 



A first upgrade of the DAMA/LIBRA set-up was performed in September 2008. One 
detector was recovered by replacing a broken PMT and a new optimization of some 
PMTs and HVs was done. The transient digitizers were replaced with new ones, having 
better performances and a new DAQ with optical read-out was installed; since October 
2008 DAMA/LIBRA is again in operation. Data of two further annual cycles are at hand. 

Considering the relevance to lower the software energy threshold of the experiment, in 
order to improve the performance and the sensitivity of the experiment and to allow also 
deeper corollary information on the nature of the DM candidate particle(s) and on the 
various related astrophysical, nuclear and particle Physics scenarios, the replacement of 
all the PMTs with new ones with higher quantum efficiency has been planned and work 
is in progress. DAMA/LIBRA will also study several other rare processes as done by the 
former DAMA/Nal apparatus in the past [12] and by itself so far lfl8ll . 



REFERENCES 



1 . R. Bernabei et al., // Nuovo Cim. A112, 545-576 (1999). 

2. R. Bernabei et al., Phys. Lett. B389, 757-766 (1996); R. Bernabei et al., Phys. Lett. B424, 195-201 
(1998); R. Bernabei et al., Phys. Lett. B450, 448-455 (1999); P. Belli et al., Phys. Rev. D61, 023512- 
9 (2000); R. Bernabei et al., Phys. Lett. B480, 23-31 (2000); R. Bernabei et al., Phys. Lett. B509, 
197-203 (2001); R. Bernabei et al., Eur. Phys. J. C23, 61-64 (2002); P. Belli et al., Phys. Rev. D66, 
043503-17 (2002). 



3. R. Bernabei et al., Eur. Phys. J. C18, 283-292 (2000). 

4. R. Bernabei el al., La Rivista del Nuovo Cimento 26 n.l, 1-73 (2003). 

5. R. Bernabei et al., Int. J. Mod. Phys. D13, 2127-2159 (2004). 

6. R. Bernabei et al., Int. J. Mod. Phys. All, 1445-1469 (2006). 

7. R. Bernabei et al., Eur. Phys. J. C47, 263-271(2006). 

8. R. Bernabei et al., Int. J. Mod. Phys. All, 3 155-3 168 (2007). 

9. R. Bernabei et al., Eur. Phys. J. C53, 205-213 (2008). 

10. R. Bernabei et al., Phys. Rev. D77, 023506-9 (2008). 

11. R. Bernabei et al., Mod. Phys. Lett. A23, 2125-2140 (2008). 

12. R. Bernabei et al., Phys. Lett. B408, 439-444 (1997); P. Belli et al., Phys. Lett. B460, 236-241 
(1999); R. Bernabei et al., Phys. Rev. Lett. 83, 4918-4921 (1999); P. Belli et al., Phys. Rev. C60, 
065501-7 (1999); R. Bernabei et al., II Nuovo Cimento A112, 1541-1552 (1999); R. Bernabei et al., 
Phys. Lett. B515, 6-12 (2001); F. Cappella et al., Eur. Phys. J.-direct C14, 1-6 (2002); R. Bernabei et 
al., Eur. Phys. J. A23, 7-10 (2005); R. Bernabei et al., Eur. Phys. J. A24, 51-56 (2005); R. Bernabei 
et al., Astrop. Phys. 4, 45-54 (1995); R. Bernabei, in The identification of Dark Matter, World Sc. 
Pub., Singapore, 1997, pp. 574. 

13. P. Belli et al., Astropart. Phys. 5, 217-219 (1996); P. Belli et al., Nuovo dm. C19, 537-544 (1996); 
P. Belli et al., Phys. Lett. B387, 222-226 (1996); Phys. Lett. B389, (1996) 783 (err.); R. Bernabei 
et al., Phys. Lett. B436, 379-388 (1998); P. Belli et al., Phys. Lett. B465, 315-322 (1999); P. Belli 
et al., Phys. Rev. D61, 117301-4 (2000); R. Bernabei et al., New J. of Phys. 2, 15.1-15.7 (2000); 
R. Bernabei et al., Phys. Lett. B493, 12-18 (2000); R. Bernabei et al., Nucl. Instr. & Meth A482, 
728-743 (2002); R. Bernabei et al., Eur. Phys. J. direct Cll, 1-8 (2001); R. Bernabei et al., Phys. 
Lett. B527, 182-186 (2002); R. Bernabei et al., Phys. Lett. B546, 23-28 (2002). R. Bernabei et al., in 
the volume Beyond the Desert 2003, Springer, Berlin, 2003, pp. 365; R. Bernabei et al., Eur. Phys. 
J. All, sOl 35-41 (2006). 

14. R. Bernabei et al., Astropart. Phys. 7, 73-76 (1997); R. Bernabei et al., Nuovo dm. A110, 189-195 
(1997); P. Belli et al., Astropart. Phys. 10, 1 15-120 (1999); P. Belli et al., Nucl. Phys. B563, 97-106 
(1999); R. Bernabei et al., Nucl. Phys. A705, 29-39 (2002); P. Belli et al., Nucl. Instr. & Meth A498, 
352-361 (2003); R. Cerulli et al., Nucl. Instr. & Meth A525, 535-543 (2004); R. Bernabei et al., 
Nucl. Instr. & Meth A555, 270-281(2005); R. Bernabei et al., Ukr. J. Phys. 51 1037-1043 (2006); 
P. Belli et al., Nucl. Phys. A789, 15-29 (2007); P. Belli et al., Phys. Rev. C76, 064603-10 (2007); 
P. Belli et al., Phys. Lett. B658, 193-197 (2008); P. Belli et al., Eur. Phys. J. A36, 167-170 (2008); 
P. Belli et al., Nucl. Phys. A826, 256-273 (2009). 

15. P. Belli et al., Nucl. Instr. & Meth. A572, 734-738 (2007); P. Belli et al., Nucl. Phys. A806, 388-397 
(2008); P. Belli et al., Nucl. Phys. A824, 101-114 (2009); P. Belli et al., in the Proceed, of the Int. 
Conf. NPAE 2008, ed. INR-Kiev, Kiev, 2009, pp. 473. 

16. R. Bernabei et al., Nucl. Instr. & Meth. A592, 297-315 (2008). 

17. R. Bernabei et al., Eur. Phys. J. C56, 333-355 (2008). 

18. R. Bernabei et al., Eur. Phys. J. C62, 327-332 (2009). 

19. K. A. Drukier et al., Phys. Rev. D33, 3495-3508 (1986); K. Freese et al., Phys. Rev. D37, 3388-3405 
(1988). 

20. D. Smith and N. Weiner, Phys. Rev. D64, 043502-1 1(2001); D. Tucker-Smith and N. Weiner, Phys. 
Rev. D72, 063509-8 (2005). 

21. K. Freese et al., Phys. Rev. D71, 043516-15 (2005); Phys. Rev. Lett. 92, 1 1 1301-4 (2004). 

22. F. S. Ling, P. Sikivie and S. Wick, Phys. Rev. D70, 123503-19 (2004). 

23. M. Ambrosio et al., Astropart. Phys. 7, 109-124 (1997). 

24. M. Selvi on behalf of the LVD coll., Proceedings of The 31st International Cosmic Ray Conference 
(ICRC2009) Lodz, Poland, 2009, in press. 

25. see e.g. B. Sadoulet, talk at TeV Particle Astrophysics; C. Galbiati talk at WIN09. 

26. F. Nozzoli, talk at TAUP 2009 Conference, Rome, July 2009. 

27. J.H. Jenkins et al., Astropart. Phys. 32, 42 (2009). 

28. PS. Cooper, Astropart. Phys. 31, 267 (2009). 

29. R. Bernabei et al., Liquid Noble gases for Dark Matter searches: a synoptic survey, Exorma Ed., 
Roma, ISBN 978-88-95688-12-1, 2009, pp. 1-53 darXiv:0806.00TlV 2). 

30. A. Benoit et al., Phys. Lett. B637, 156-160 (2006). 



